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Abstract 

In this paper we investigate dynamics of the modified loop quantum cosmology models using 
dynamical systems methods. Modifications considered come from the choice of the different field 
strength operator F and result in different forms of the effective Hamiltonian. Such an ambiguity 
of the choice of this expression from some class of functions is allowed in the framework of loop 
quantization. Our main goal is to show how such modifications can infiuence the bouncing universe 
scenario in the loop quantum cosmology. In effective models considered we classify all evolutional 
paths for all admissible initial conditions. The dynamics is reduced to the form of a dynami- 
cal system of the Newtonian type on a 2-dimensional phase plane. These models are equivalent 
dynamically to the FRW models with the decaying effective cosmological term parameterized by 
the canonical variable p (or by the scale factor a). We demonstrate that the evolutional scenario 
depends on the geometrical constant parameter A as well as the model parameter n. We find that 
for the positive cosmological constant there is a class of oscillating models without the initial and 
final singularities. The new phenomenon is the appearance of curvature singularities for the finite 
values of the scale factor, but we find that for the positive cosmological constant these singular- 
ities can be avoided. The values of the parameter n and the cosmological constant differentiate 
asymptotic states of the evolution. For the positive cosmological constant the evolution begins at 
the asymptotic state in the past represented by the de Sitter contracting (deS_) spacetime or the 
static Einstein universe H = or H = —oo state and reaches the de Sitter expanding state (deS+), 
the state H = or H = -|-oo state. In the case of the negative cosmological constant we obtain 
the past and future asymptotic states as the Einstein static universes. 
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I. INTRODUCTION 

According to the General Relativity the initial state of the universe was a causal singular- 
ity. In this point all time-like geodesies were cut off and has no extension into the past. This 
peculiar behaviour corresponds to the blow up of the curvature IZ ^> oo and energy density 
p ^ oo in the standard Big-Bang scenario. It is however believed that singularities of this 
kind are not real states but rather artifacts resulting from the fact that classical theory broke 
down. In the vicinity of the initial singularity, when the universe reaches the Planck scale 
energy, the quantum gravity effects were proposed as a mechanism to remove the singular 
behaviour. Such a possibility was studied in the context of the String Theory leading to 
the so called pre-Big Bang cosmological scenario [l[. Presently this issue is also studied 
extensively in the Loop Quantum Cosmology (LQC)|2] where the avoidance of the initial 
singularity appeared in a natural way [3| . This approach possesses other advantages like the 
adjustment of the proper initial conditions for the inflationary phase. The avoidance of the 
initial singularity appeared here in the form of the non-singular bouncing universe y, Isj. 
Namely, the universe is initially in the contracting phase, then reaches a minimal nonzero 
volume and due to quantum repulsion evolves toward an expanding phase. It should be 
stressed that this result has been strictly demonstrated so far only for isotropic models with 
a free, massless scalar field |6|. The phenomenological dynamics of this scenario has been 
extensively studied [7|, l8|. However ambiguities in the quantization of the classical formula 
lead to some deviations. In the frame of the standard quantization scheme these ambiguities 
lead generally to the models with qualitatively similar bouncing behaviour. However like 



it was shown in 



for some special setup in the effective fiat FRW model the bounce is 



replaced by non-singular oscillations. More drastic consequences appear when we choose the 



Is 



different quantization framework. This issue was studied in [9| and it was found the existence 
of the finite scale factor curvature singularities in the co-existence with a quantum bounce. 
In this paper we continue these studies in the nonstandard quantization frameworks. We 
investigate the effects of resulting ambiguities for the dynamics of the effective FRW models 
with the cosmological constant. We present simple methods of reducing dynamical equation 
governing evolution of the models with quantum corrections to the form of a 2-dimensional 
autonomous dynamical system of the Newtonian type. We explore a mechanical analogy in 
investigation of different evolutionary scenarios. 



The organization of the text is the following. In section [TTl we investigate possible kinds of 
ambiguities in the loop quantization, concentrating on the expression for the field strength 
operator F. Next, in section IIIII we investigate resulting cosniological models with use of 
the qualitative methods analysis of dynamics. In section |V] we summarize the results. 

II. FIELD STRENGTH OPERATOR F AND QUANTIZATION AMBIGUITIES 

This section aims to summarize some key features of loop quantization of cosmological 
models and discusses arising ambiguities. In particular we extend study of the choice of 
the field strength operator F began in [9|. We restrict our consideration to the fiat FRW 
cosmological models which are given by the metric 

ds"^ = -N'^{x)dt^ + qabdx^dx^ (1) 

where N[x) is the lapse function and the spatial part of the metric is expressed as 

g„, = 5,,u^W, = a\trq^, = a2(t)<5,/a;r^^ (2) 

In this expression °g^j, is the fiducial metric and "cj^ are the co-triads dual to the triads °e", 
°uj^{°ej) = 6j where °uj^ = °ujldx"- and °ej = "e^da- Based on the triads we can construct the 
Ashtekar variables 

A^ ^ r^ + 7i^: = cV-'/' "u;l, (3) 

Ft ^ ^Ad^\e1 = pVo-'^'V^ X (4) 

where 

\P\ = aX^\ (5) 

c = ^dV,'^\ (6) 

The volume Vq is the so called fiducial cell and is introduced to perform the integration 
on the non-compact spaces. We expect that the physical results should not depend on its 
value. However, when we introduce quantum modifications in general they can depend on 
the choice of Vq. 

In loop quantum cosmology we consider effects of the inverse volume corrections and 
holonomies. In the phenomenological considerations of the quantized FRW k = model it 



is appropriate to take into account quantum modifications arising only from tlie field strength 
expressed in terms of holonomies. In some special case {n = —1/2) this kind of modification 
does not depend on the choice of the fiducial cell Vq. This simplifies interpretation of results 
and make the correspondence with classical physics easier. This issue Js related to the 
question of the loop quantization uniqueness of the cosmological models |101]- Because there 
is no defined length scale in the fiat FRW models the Vq dependence seems to be unphysical. 
This problem is not present in the FRW fc = ±1 models, where the length scale is tied with 
the intrinsic curvature ll[ |. 



The loop quantization requires to replace the classical variables by the quantum operators 
expressed in terms of the holonomies and fluxes. In the FRW models holonomies reduce 
to almost periodic functions Ylj^j^^^^^^"^ ^^^ fluxes reduce to variable p multiplied by the 
kinematic factor 12j. Then algebra of variables p and e*'^'^/^ is quantized. The quantum 
Poisson bracket of elementary operators is 



1/J.C 

j9, e^" 



87r/p,7'^ 



(7) 



and these operators act in the Hilbert space ^^(IRBohr, dfiBohr)- The basis of this space can 
be chosen as a set of eigenstates of the operator p 



p\fi) = fi' 
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PI 



1/^) 



which fulfils the normalization condition (/^|/^') = S^^^r. The action of the operator A^ 



is following 



6 2 1^ ) = l/i +/i). 



(8) 



e 2 



(9) 



We have defined a quantum kinematic sector of the loop quantum cosmology. Now with 
use of these rules we aim to quantize the scalar constraint 



n 



G 



d^xN{x) 



E-E] 



where the extrinsic curvature is defined as 

1 



e^hF',,-2{l + ^^)KlK, 



(10) 



K. 



ah 



2N 



[qab - 2D(^aNb)\ 



(11) 



which corresponds to if* := KabC^. The scalar constraint quantum operator should be 
expressed in terms of the elementary operators p and N = e^~ . In our considerations we 



restrict ourselves to quantization of the field strength F^^. The rest of the terms lead to 
quantum operators expressed in term of the inverse volume operator. Detailed analysis of 
the scalar constraint quantization can be found in [l3|. 

As we have said our main interest is to find a formula for the field strength operator F. 
The standard approach to introduce this operator is to consider a closed path a = Dj^ G S 
and holonomy along this loop. Such a holonomy can be written as 






(m)l(m)l(m)-1l(m)-1 



hrhf'h 



h\ 



(12) 



where particular holonomies are 

/,(/^) ^ gr.^c ^ J ^Qg j^/f^ j ^ 2r, sin (^) . 



(13) 



With use of the above expressions one can derive formula 



tr 



ab 



-2- 



^(/^[£ 



^^^v^/' 






g(/ 



(14) 



Based on this expression we can define the quantum field strength operator. However, 
we should discuss some related issues. Firstly, quantization of the area operator in the Loop 
Quantum Gravity leads to the existence of the area gap Armin = A 14| . This means that the 
limit Ar^. . — > of the area enclosed by the loop does not exist. We should stop shrinking 
the loop at a minimal area which is related to limiting yU -^ /i value. This process is however 
ambiguous and generally limiting [i can have a form 



-^ = iW 



(15) 



where n E [—1/2,0) [15[ and ^ is dimension-full and positive real parameter. In particular, 
from this expression we can recover known /2-scheme for which 



/^ 




(16) 



In this case the area gap A corresponds to the physical area of the loop |p|/i^. The range 
of the n parameter written above can be obtained from the considerations of the lattice 



states IGj. Further restrictions for value of this parameter come from considerations of 



the anomaly cancellation in the inhomogeneous cosmo^ 



graviton effective mass, giving —0.1319 > n > —5/2 17l |. 



ogical models and positivity of the 



The properties considered should be now taken into account during the construction of 
the operator F . The standard procedure is to cut-off the terms 0{^'^) and perform the hmit 



Ar — > A in expression (IT^ . what gives [18|] 



tr 



r, (/.[^ - I 



F^, = -2 hm „ ,„ 



This is the well defined operator on the Hilbert space L^(MBohr) c^/^Bohr)- Namely, due to 
relation sin {jic) = ^ {e^P"^ — e"*'"'^) it can be expressed in terms of the operator N and due 
to the power-low p,{p) dependence this part can be expressed in terms of the operator p. 
The proper classical behaviour is also obtained, it means that no information about classical 
expression is contained in the terms (9(/i^). These terms simply vanish in the classical limit 
fi —>■ 0, 0(/i^) —>■ 0. Another motivation to consider (ITTll as a proper quantum operator is 
the resulting dynamics. Namely, this quantum operator leads to the avoidance of the initial 
singularity. 

However we can ask what happens when a part of the 0(/i^) contribute to the definition 
of the operator F7 The resulting operator has proper classical behaviour, because the 
contribution from O(yU^) vanishes in the classical limit. However on the quantum level 
the operator F contain some new terms which influence quantum dynamics. It is worth 
to stress this property. This additional modification term, which influence dynamics in 
the quantum regime, does not correspond to any contribution on the classical level. The 
possibility of existence of such terms lead to the ambiguity in the definition of the operator 
F. However not all modifications which fulfil well the classical limit are allowed. Another 
property which must hold is that the resulting operator must be defined on the Hilbert space 

In the paper [9| the expression for the F operator including the term 0{fi'^) was calculated 



1 — -i/l — |sm (uc) 

fpk _ k V ■^ o «o J /-|Q\ 



10V'" 



r»/o 



Expanding the square we obtain 

9 /- \ 

sin (/ic) 



1 — -1 / 1 — o sin {jj,c) 



3^ 



1 sin"^ (/ic) 



+ 0— 1^ + ---- (19) 



We see that the first term in this series corresponds to the standard field strength operator. 
Other parts vanish in the limit /i — i> and have no classical meaning. With use of the 
formula sin (/ic) = ^ (e*^'^ — e"*'^'^) we can express the series in terms of the operator N . 



In the next section we apply the modified expression for the field strength operators 
investigated above in the cosmological models. Our goal is to find whether such modifications 
lead to singular behaviour and how the bouncing universe scenario is modified. We will 
restrict to the fiat FRW models with the free scalar field and the cosmological constant 
content. Our analysis concentrates on the investigation of the effective equations rather 
than strict quantum dynamics. So we obtain information only about evolutions of the mean 
values of the canonical parameters. The information about dispersions is unavailable in such 
an approach. We will consider both full modifications from the term 0{^'^) as well quadratic 
modifications from the expression (TT9|) . 

III. QUALITATIVE ANALYSIS OF DYNAMICS 

The dynamical equations resulting from loop quantum cosmology corrections cannot be 
integrated in the exact form and that is why the qualitative methods of differential equations 
can provide us with a powerful tool for investigations of evolutional paths of the models under 
considerations. 

The main advantage of using these methods lies in possibility of examine of all solutions 
for all admissible initial conditions in the geometric language of a phase space (space of all 
states at any time). As a result we obtain description of global dynamics of these models, 
i.e., the phase portrait from which we obtain their clear stability analysis. In the dynamical 
systems framework the crucial role plays the critical points (for which right-hand sides of 
dynamical system vanish) and the trajectories joining them. From the physical point of view 
the critical points represent stationary states of the system or asymptotic states. On the 
other hand trajectories represent evolution of the system. The famous Hartman-Grobman 
theorem |l9J tell us that linearization of the system at this critical points (non-degenerate) 
contain all necessary information about its behaviour in neighbourhood of this critical point. 
From the linearization matrix we obtain all necessary information concerning the type of 
the critical point as well as its stability. 

The dynamics of cosmological models under considerations can be reduced to the 2- 
dimensional dynamical systems of the Newtonian type with the total energy E and the 
potential V. As it is well known [19| for such systems only two types of critical points 
are admissible, namely a saddle (if eigenvalues of the linearization matrix are real of oppo- 



8 



site signs) or a centre (if eigenvalues of the linearization matrix are purely imaginary and 
conjugate). While the first type of critical points is structurally stable, the second one is 
structurally unstable, i.e., small changes of its right-hand sides disturb the global phase 
portrait which is defined modulo a homeomorphism which establishes the equivalence of 
the global dynamics. Of course the motion of the system in the configuration space takes 
place in the domain admissible for motion determined by condition of non-negativity kinetic 
energy form, E — V is non-negative. The critical points of the system are always located 
at the extrema of a diagram of the potential function (if we have an inflection point on the 
diagram then the corresponding point is degenerated). If in the diagram of the potential 
function we have local maxima then they correspond to saddles on the phase portrait. In 
turn if we have local minima on the diagram then in the phase space the centre type of the 
critical point is present. Note that for systems under consideration all critical points in the 
phase space must lie on the p axis and they represent static Einstein universes. Apart from 
the critical points at the finite domain there are present physically interesting critical points 
at infinity. Therefore for the full study of djTiamics of the system the analysis of behaviour 
of trajectories at infinity is necessary. For this aims we complete the plane by adding to 
it a circle at infinity via the Poincare sphere construction or equivalently by introducing 
the projective coordinates (the construction of the projective plane) (l9|. In our case on 
the circle at infinity critical points are located which represents the de Sitter universes if 
n = -1/2. 

Due to the particle-like description of dynamics in the notion of classical mechanics we 
obtain a simple description of evolution in terms of particle moving in the potential well. 
We can classify all evolutionary pats in both configuration and phase spaces. It is important 
that the dynamics is fully characterized by a single potential function V{p). If we have 
the global phase portraits of the cosmological models then one can discuss how different 
models with a desired property, say a singularity, bounce, acceleration are distributed in the 
ensemble. Are they typical or generic? The answer for these questions are interesting in 
cosmological context when problem of initial conditions is crucial. 



A. (^(/U^) modifications 



In the case of 0(^^) modifications tlie phenomenological Hamiltonian takes a form 

9 



H 



phen 



2^7^/2 



2r,2 



4 



1- Wl--sin2(/ic) ]^M+\P\ 



3/2 



P 



(20) 



where matter part contains the free scalar field and the cosmological constant 



pI 

2W 



A 

+ -. 

K 



(21) 



The Hamiltonian [2D] traces the effects of the 0(/i^) modifications. However it must be 
emphasized that it is not necessarily the true effective Hamiltonian. It mean that solutions 
of the Hamilton equation / = {/, Tiphen} not necessarily traces the mean values of the 
quantum operator (/). This issue must be verified by the full quantum treatment. 
From the Hamiltonian equation we obtain evolution of the canonical variable p 

2vbl sin {j2c) cos {j2c) 



P = {P, ^phen} 



7/x 



(22) 



— I sin (fic) 

Based on this result and the Hamiltonian constraint Tiphen = one derives the modified 
Friedmann equation in the form 



H' = 



It. 
4p2 



3 V 3pc 



3 1 

4 + 4 



3 pc 



where 



Pc 



Now we make the following change of the variable x 
the form 



(23) 



(24) 



P 



l-n 



and rewrite equation fl23l) to 



X 



3 Pc 



2 (l-n)2Kp^^2 /^ 



3 Pc 



-xHi 



1 A l + 2n . 

3Ao 



1 a 
'2T0 



(25) 



where a 



K^2 



and Ao 



iW 



. Due to the symmetry of the change p — *■ —p we can consider 



only positive branch of p, as it was done above. In such a approach negative solutions 
can by simply obtained multiplying positive solutions by minus one. Performing the time 
transformation 



d 
d^ 



\ 



1 



I 



IP. 

3 Pc 



d 



(1 — ny Kp -. 



3 1 



3 Pc 



dt 



(26) 
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and introducing the potential function is in the form 



and define 



l/(x) = -\.Hl - —-**') (27) 



'^-li (-) 



-— ( a + '-x^^ 1 =x^ -x'"^ > 0. 



we rewrite the modified Friedmann equation in the point particle form 

-x'^ + V{x) = E. (29) 

In Figs [1], [21 [3] and m we present the phase portraits for various values of the cosmological 
constant A and the parameter n for the model with modification 0{fi^). 

In ( l26l) the term under square root needs to be positive. In our new variables {x, x') the 
condition p > fl2Tl) corresponds to 

If A ^ 

The physical region in the phase space is defined by the condition x'"^ < x"^. In the phase 
portraits this region is labelled as white. The time transformation (126|) is singular at 1 — |— = 
which, in the new variables, corresponds to x'"^ = |x^. Location of such points in the phase 
plane denote curvature singularities of the original system (iJ^ = cxd see ( l23l) ) and are labelled 
as a solid blue line (electronic version). 



B. Quadratic modifications from ©(/x^) 

The employ of the quadratic contribution from the 0{fi'^) studied in the section HTl lead 
to the phenomenological Hamiltonian in the form 

Hphen = |r^ ( sin^ (fic) + ^ sin^ (fie)) + \pf^p (30) 

where matter content is the same as in the previous case and is given by ( 12T|) . In the same 
way like in the previous case we aim to calculate the effective Friedmann equation and then 
rewrite it in the point-like particle form. To achieve it we derive 



p = {p, Hphen} = — sin {nc) cos (/ic) (l + ^ sin^ (pc)) (31) 

7/i o 
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FIG. 1: The phase portrait for A = and arbitrary n in the case A type modifications. The 
shaded region is physicahy forbidden. The dark gray hne (blue el.v.) denotes phase space location 
of H^ = CO. The light gray line (red el.v.) which is a border of the physically allowed region, 
denotes phase space location of H^ = 0. The zero energy level trajectories joining critical point at 
the finite domain with critical points at infinity divides the phase space in to two region occupied 
by the scalar field (a > 0) and the phantom scalar field (a < 0) trajectories. Only trajectories with 
a > are physical and for all admissible initial conditions they tend to the final state of H^ = 
for x = oo (p = cxDora = cxd). 



and with use of the Hamiltonian constraint we obtain the modified Friedmann equation 

5 



3 V 3pc 



3^A + li 



1 + ^/1 + ti; 



Changing of variable x = p^ " we rewrite equation 

1 



to the form 



1 + JI + 



4 p 



x" i 3 -^ "^ V -^ "^ 3 pc 12/^ 4 A i+2u 



- — ^ j \ -i- ^ 1 — n 

2 (l-n)2fi:p5_3 /i + Ij^ 2 V 3Ao 



3 pc 



Performing the time transformation 



d 
dr 



^ 3 Pc d 



^(i--)^-P5-3yrT|^ 



dt 



2f- 
Ao 



(32) 



(33) 



(34) 
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FIG. 2: The phase portraits for n = —1/2 in the case A type modifications, a) For < — < 2 the 
scalar field (a > 0) and the phantom scalar field (a < 0) trajectories are physically allowed, but 
their evolution is different; the a > trajectories start from anti de Sitter state and go through 
hyper-deflation state then bounce appears and next we have hyper-inflation state, the final state 
is de Sitter; the phantom trajectories for x' > start from H^ = and tend to de Sitter state; b) 
For I < 7p < 3 the situation is quite similar but the state of hyper-inflation (or hyper-deflation) 
appears only for phantom trajectories. 
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FIG. 3: The phase portraits for n = —1/2 in the case A type modifications: a) For — = 2 ^^^ 
trajectories with a > after bounce tend to the critical state of H^ = cxd for x = cxd, the phantom 
(a < 0) evolution for x' > starts at H^ = state and tends to the same final state of H^ = 00; 
b) For A < only the scalar field trajectories are physically allowed (a > 0), their evolution starts 
at H^ = and after hyper-deflation, bounce and hyper-inflation state ends at H^ = 0. 
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FIG. 4: Phase portraits for a) A > and n = 0; b) A < and n = 0; c) A > and n = -1/4; d) 
A < and n = —1/4 in the case A type modifications. Note that the cases n = and n = —1/4 
are topological equivalent. 



we obtain Hamiltonian of the point-like particle of the unit mass 

-x'^ + V{x) = E 
where the potential function is given by 

T^/ N 1 2/. 4 A l + 2n\ 



(35) 



(36) 
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FIG. 5: The phase portrait for A = 0, arbitrary n and the type B quadratic modifications. The 
shaded regions are physically not allowed. The borders of these regions denote the phase space 
location of the state H^ = 0, only the scalar field (a > 0) trajectories are located in the physical 
domain. For x' > they start from the state H = and tend to the final state of -ff = at a; = co. 
The states H = represent the static Einstein universes which are reached asymptotically. 



and total energy 



Ao 



(37) 



Figs [SI [6], [3 and [8] represent the phase portraits for model with this type modification. 

For this type of corrections the condition p > corresponds to x''^ > x^ . There is another 
condition for positivity of the Hubble function (152]) and the time transformation fl34p which 
is I > \ h- + |- . In the new variables (x, x') this is x'"^ < ^x"^. The physical domain in the 
phase space for the type B corrections is defined by 

x'^ < x'^ < — x^. 
- 9 

This region is labelled as white on the phase portraits. 



16 




1/2 and the type B quadratic modifications: a) For 7^ = 3 



FIG. 6: The phase portraits for n - 

only a < trajectories are physical and for x' > they start form the H^ 

and go to H^ = for x = cxd; b) For < 

start form H^ = with finite x and tend to the de Sitter state at infinity. 



Ao 
state for finite x 

^ < I both types of trajectories are allowed and they 
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FIG. 7: The phase portraits for n 

Ao 



-1/2 and the type B quadratic modifications from 0{fj, ): a) 



^>l>b)-|<A<o. 



C. Exact solutions for the type A and B modifications 

It is interesting that for both types of modifications we can find exact solutions for x{t) 
or p{t) for n = —1/2 and n = 0. In the case of n = —1/2, a > the modified Friedmann 
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FIG. 8: The phase portraits for model with the type B quadratic modifications for n = 0: a) and 
b), n = —1/4: c) and d) and A > 0: a) and c) (left panel), A < 0: b) and d) (right panel). The 
physical trajectories lie in the domain labelled as white. 

equation for the type A modification takes the form 



/c^y ^ _^ 2 / _ 1 A_ 

\drj Ao V 3Ao 



Tlie pliysical solutions of this equation correspond to the case 

A 

Ao 



(38) 



(39) 
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Performing the direct integration of the equation (138|) we obtain solution 



XlT 



3a 



3An-A 



cosh 



1 A , 



or 



p(r) 



3a 



1/3 



COS 



h2/3 



1 A , 
^-3A;(^-^°^ 



(40) 



(41) 



,3Ao-A, 
for the canonical variable. 

In the case of n = —1/2, a < (phantom) physical solutions are permitted for A G M. 
For A/Aq < 3 we obtain 



x{t) 
for A/Aq = 3 we obtain 



—3a 
3Ao-A 



cosh 



1 A , 



x[t} =Xo+ J-—T, 



a 



for A/Aq > 3 we obtain 



x(r) 



—3a 

A-3An 



sm 



1^ U ^ 



(42) 



(43) 



(44) 



Solutions for p(r) are obtained taking p(r) = x^/^(r). 

For n = the modified Friedmann equation takes the form 

(dx\ a 2 1^3 

The solution of this equation is expressed in terms of the Weierstrass p-function 

Ao 



x{t) = ^ [1 - 12p(r - To; ^2, ^3)] 



(45) 



(46) 



where 



92 - -, 



^3 



1 / A\ a 



(47) 
(48) 



216 144 \AoJ Ao 

For the A < 0, a > this solution represents a bouncing type of evolution. In the case 
A > 0, a > physical solutions are present when additional relation 



4 a /A 

3 ^ a; Ia;^ 



(49) 
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TABLE I: The Hamiltonian and derivative of the potential at the critical point for both types of 
corrections 



Hamiltonian 
potential 



critical points : 
arbitrary n 

-1/2 < n< 



-Ha 

V{x) 



lr'2 



+ Vix) 



1^2 



iA 

3Ao 



_l_a_ 
2Ao 

l + 2n 

X 1-" 



xo = 0, 






1 



l + 2n 
, 1-1 



Xo 






2(l-n)fi,A>0, 

_ l+2n 



Zfl_ 



1-n 



B 



Wb 



y(x) 



-x'"^ 



1^2 



+ y(x) = 2^ 



1 + 



4_A_ 
3Ao 



X 1-" 



Xo = 0, 
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is fulfilled and represents an oscillating type of evolution. In this case solutions for x{t) and 
pij) are equal, p(r) = x(r). 

Analogous solutions for the type B modifications can be also obtained with the help of 
transformation A -^ — 4A and a -^ —4a (see Table H]). 

IV. GENERAL PROPERTIES OF THE PHASE PORTRAITS 

The dynamical behaviour of the models considered is fully determined by the single 
potential function fl271) or fl36|) . Its diagram characterizes a type and location of the critical 
points of the dynamical system 



X 



y 



" dx 



(50) 



System (150|) is of the Newtonian type and at the finite domain of the phase space there can 
be only a saddle or a centre type critical points (or degenerated which are not interesting 
because we are only looking for generic cases). The type of a critical point is determined 
from the characteristic equation of the linearization matrix, namely 



y 



dx"^ 



(51) 



xo 



calculated at the critical point. If 



dx^ I 



xo 



> eigenvalues are real of opposite signs and we 
d^v\ 



have a saddle type critical point, and in case if — ^^ | < eigenvalues are purely imaginary 
of opposite signs and a critical point is of a centre type. 

For the model with type A modifications (see also Table [T]) 



21 



n = -1/2, xo = 0: A = ±^l - i^ 
a saddle for ^ < 3 
a centre for x- > 3 

Ao 

-1/2 < r2< 0: 

xo = : A = ±1 a saddle 

Xn'"" = 2(1 - n)4a and A > 0: A = ±i^l^ a centre 

For the model with type B modifications (see also Table [T]) 



. n = -l/2, Xo = 0: A = ±^1 + |^ 
a saddle for x- > — f 

Ao 4 

a centre for ^ < — | 

• -1/2 <n< 0: 

Xo = : A = ±1 a saddle 

l + 2n I 

Xq-" = -i(l - n)^ and A < 0: A = ±iJ^ a centre 

Therefore all information about the critical points and their type contains the single 
potential function. 

Figs [IHl] show representative cases of the global dynamics of the system of equations fl29|) 
and (l27|l . Figs [5H8] show representative cases of the global dynamics for type B corrections. 

We obtain different phase portraits for basic model parameters (n, A), where —1/2 < n < 
and the cosmological constant belongs to the some characteristic intervals. For comparison 
we also present phase portraits for the case of the vanishing cosmological constant Fig. [T] for 
the case A and Fig. [5] for the B case. In the case A we can observe only critical point of a 
centre type located at the x axis and critical point of a saddle type at the origin of the phase 
space (in principle at the origin a centre type critical point for A/Aq > 3 can be present, but 
this case in not interesting because of it's structural instability). Along the trajectories it is 
measured the new time variable which is a monotonic function of the original cosmological 
time but this flow can be incomplete, i.e., inflnite cosmological time is mapped on the finite 
interval measured along the trajectory. 
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In the case of the positive cosinological constant and the case A modifications two char- 
acteristic evolutional scenarios are realized: 

• A trajectory starts from the contracting deS_ stage and finishes at the expanding deS+ 
stage. In the middle state the trajectory is undergoing the bounce in the vicinity of 
the intersection with the p-axis. In other words for such a evolution deS_ is the initial 
state and deS+ is the final state (see Fig. [9l). 

• There is new non expected behaviour for the positive cosmological constant - emerging 
of oscillatory behaviour - without initial and final singularities and, in special case, 
without hyper-infiation and hyper-defiation states (see Fig. (TD]) for 



a 



1^ 2 ~ 1/. , ^ N //. ^ ^ \2t+2S 



-pi>-{l + 2n)Ul-n)- 



6 



6 



Ar 




hyper-dcflation 



classical evolution 



classical 
singularity 



FIG. 9: The schematic picture of the evolution with the positive A for the monotonic models. The 
green curve represents the energy density. It is worth to note that this energy density is finite 
during whole evolution, even during transitions through singularities. The blue curve represents 
the schematic evolution of the scale factor for the investigated model. The dashed curve represents 
the classical evolution which is not realized in the presented quantum model. 

What it is interesting in these scenario is the universe is undergoing a phase of hyper- 
infiation in the vicinity of a singularity. This singularity is a curvature singularity at which 
the scale factor (as well as density) is finite but the Hubble function assumes infinite value. 
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• - finite scale factor singularities 



recollapse 




time 



FIG. 10: Generic evolutional scenarios for bouncing-oscillatory models with A > and the value 



of critical energy constant is Ec = — ]^(1 + 2n)l (1 — n)— 1 " : a) the oscillating model without 
the singularities {E < Ec); b) the oscillating model with double singularities at finite scale factor 
per period {E = Ec); c) the oscillating model with four singularities at the finite scale factor per 
period {E > Ec). 



We mark a blue line in the phase portraits which shows the singularities reached by different 
trajectories. The algebraic equation for this line is independent of the value of the cosnio- 
logical constant. But note that there are trajectories which do not intersect this line at all, 
they intersect it only at one point and at least they intersect it at two points. These three 
cases have been illustrated schematically in Fig. [TOl li E < Ec singularities are avoided. 

The new evolutional scenario offered by the quantum correction is oscillating behaviour 
in universes with double singularities during the expansion (and contracting) period for the 
positive A. For other propositions of oscillatory models with bounce and phantom see 20|. 

All emerging evolutional scenarios for the origin of the universe from the bounce are in 
principle testable because we have the Hubble function in respect of the redshift, which in 
order can be starting point for cosmography. In principle if the size of the universe is going 
to infinity there are three types of exits: H is going to infinity, H is approaching to the de 
Sitter state or H is going to the zero. All these long term behaviour one can find in Fig. [9l 

The quantum corrections considered are dynamically equivalent to the effects of the 
decaying cosmological constant term with the effective equation of state coefficient param- 
eterized by the scale factor or the redshift. This fictitious fluid which mimics dynamical 
effects of the quantum corrections can be calculated in the exact form for the both types of 
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FIG. 11: The dependence of the effective equation of state parameter on the redshift z for both 
cases A (top) and B (bottom) of quantum corrections for n = —1/2 and A > (see phase portraits 
Fig. [2^ and Fig. [6)3). In both cases energy density approaches the cosmological constant w^.^ = — 1 
case at the long time evolution. Top panel: in the type A correction the evolution begins at 
the bounce for the redshift value Zb, next the curvature singularity is present at z^ and then WeS 
approaches the —1 value. Bottom panel: in the type B correction the evolution begins at zq when 
H^ = 0, next the crossing of the phantom divide line is present at z^ (marked as a dot) and then 
again w^s approaches —1. 
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corrections. Using the general formulae 



Weff = -1 



IdlnH^ 



3 din a 

and modified Friedmann equations fl23p and (1321) we obtain following forms of the equation 
of state parameter: for the type A correction 
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(52) 



and for the type B correction 
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(53) 



(54) 



where 

Pc Ao V 3. 

Equation fl54|) gives us opportunity to calculate the redshift distance to the characteristic 
moments of the universe evolution. For the type A correction the bounce occurs when 
p/pc = 3 and using ( l54l) we have that at the scale factor value fulfilling relation 

the bounce is present. In the special case of n = —1/2, the value of the redshift at the 
bounce is 

Another characteristic moment of the model evolution is the curvature singularity appearing 
at the redshift (for n = —1/2) 



VAo 



-'•/^(^^£'- 



(56) 



For the type B corrections evolution begins when p/pc = 4/3, that is, for n = —1/2 at 

(57) 



zo 



3 a V3 



A 

a; 



There is another important moment in the evolution appearing only for the type B correction, 
namely the crossing of the phantom divide line (Fig. [11] bottom). In this case the during the 
evolution of the universe Wes starts from a large negative value, crosses the line Wes = — 1, 
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then reaches the inaxiinuin and now decreases to asymptotically to minus one. Solving 
Eq. fl53|) for w^g = — 1 we find that p/pc = |(1 + 2\^) and for n = —1/2 at the redshift 
value 



-1 + 



lAn /I.. „ /TT. A 



1(1 + 2v^)--) (58) 



3 a V8' Ao 

the crossing of the phantom divide line happens. 

V. CONCLUSIONS 

In this paper we have considered the infiuence of the modifications to the field strength 
expression for the effective dynamics in the loop quantum cosmology. We have studied 
the sensitivity of evolutional scenarios on the choice of model parameters {n,A). For the 
classification of evolutional paths in arising the loop quantum cosmology models it is natural 
to adopt the dynamical systems methods which offers full classification of dynamical outcome 
for all admissible initial conditions. We consider two types of quantum modifications labelled 
as modification (A) C(/i^) and the quadratic modifications from 0{p'^) (B) (for details see 
Section HTj) . 

From mathematical point of view the dynamics is represented by 2D dynamical system 
of the Newtonian type. From the cosmological point of view quantum corrections are man- 
ifested by decaying part of the cosmological constant A term. Main advantage of using 
dynamical systems methods lies in the possibility of studding complexity of dynamics. In 
the geometric language of the phase space we can analyze the asymptotic stability of the 
solutions. From this analysis we obtain that a generic feature of dynamics is emergence of 
the bounce instead of the initial singularity. We can find easily beginning and final states of 
dynamics in its long time behaviour using the phase space compactification with the circle 
at infinity. For the positive cosmological constant and (A) type corrections, in principle, 
there are three types of outcomes (and incomes respectively because the dynamics is invari- 
ant under the refiection symmetry p -^ —p), namely the expanding de Sitter state deS+, 
H = oo state or the Einstein static universe H = 0. 

The new feature of dynamics within models with (A) type corrections is appearance of 
singularities at the finite scale factor and the super-accelerating phase of expansion in the 
vicinity of this singularity. However, we demonstrate that this singularity can be avoided 
for some choice of the critical value of total energy Ec- 
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Another new phenomenon is the emergence of oscillating models for positive value of the 
cosmological constant with single or double curvature singularities for the finite scale factor. 
We demonstrate that the dynamics with (A) and (B) a type of corrections are dual and 
after the transformations A -^ — 4A and a -^ —Aa one can obtain the evolution of models 
with the massless scalar field from evolution of the phantom scalar field (formally this means 
that we consider a opposite sign of total energy and the cosmological constant). The main 
difference between the models with positive and negative cosmological constant is that in 
the latter case incomes and outcomes always represents the static Einstein Universe. 
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